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ABSTRACT 

We present a detailed analysis of the surface brightness distribution of the brightest cluster galaxy (BCG) 
in each of 24 galaxy clusters at 0.03 < z < 0.13. We use two-dimensional profile fitting to model the surface 
brightness out to r = 300 kpc for each BCG, comparing r'/"*, r^'", and double r'/'* models. We obtain statistically 
superior fits using a two component model consisting of a pair of r'/"^ profiles with independent scale lengths, 
ellipticities, and orientations. The two component model can simply reproduce the observed position angle and 
ellipticity gradients, which cannot generally be explained purely by triaxiality. The inner component of our 
two component model has properties similar to a typical massive elliptical galaxy and is clearly associated with 
the BCG. The outer component is 10-40 times larger in scale, has '^ 10 times the total luminosity of the inner 
component, and exhibits a steeper < /i > -r^, relation than that of the elliptical fundamental plane. We interpret 
this outer component as a population of intracluster stars tracing the cluster potential. The two components are 
strongly aligned (\Ad\ < 10°) in roughly 40% of the clusters. When they are not aligned, the components tend 
toward high degrees of misalignment, suggesting that accretion of infalling material may change the orientation 
of some BCGs for a time. The extent of the outer component and its similar elongation to published cluster 
galaxy distributions indicates that the evolution of the intracluster light is tied to the cluster as a whole rather 
than to the BCG. 
Subject headings: galaxies: clusters: general — galaxies:cD, formation, evolution, fundamental parameters 



1. INTRODUCTION 

One of the most distinctive characteristics of brightest 
cluster galaxies (BCGs) is that many exhibit shallower ex- 
tended surface brightness profiles than their lower mass 
ell iptical counterparts. This t rait w as first recognized 
by Mathews, Mor2an& Schmidt (1964), who observed that 
some BCGs emit more light at large radii than would be ex- 
pected for an r'/^ profile. This excess emission, which is of- 
ten lumped into the term "cD" envelope,^ was confirmed by 
[Oemler ( 1973, 1976) and Schombert ( 1986, 1987, 1988) in a 
subset of BCGs. 

What is the nature of this excess luminosity at large 
radii? Schombert, who considered two-component models 
when fitting surface brightness profiles to r '--^ 200- 1500 
kpc, argued that the extended light is physically dis- 
tinct from the central galaxy. In this picture, the excess 
light component is thought to arise from a population 
of intracluster stars produced by th e tidal stripping and 
disruption of cluster galaxies (c.f. 'Gallasher & Ostrikei' 
19721 iDressled 11 979; Richstone & Malamuth 1983; 
Miller 1983; Malamuth & Richstone 1984; Merritt 1984; 
Spmmer-Larsen, Romeo, & Portinari 2004; Murante et al. 
2004; Willman et al. 2004). Simulations predict that 
this intracluster light (ICL) is comprised of old stars 
that are dynamically distinct from t he sta r s in the BCG 
iSommer-L arsen. Romeo. & Portinari" "2004*, iMurante et al.l 
12004ft and that production of intracluster stars is an ongoing 
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^ This emissi on is differentiated into categories by some experts (see 

ISchomberliT9^ . 



process JWillman et al.l 120041) . If correct, this model im- 
plies that the properties of the ICL, including its fractional 
contribution to the total cluster luminosity and its spatial 
distribution, are sensitive probes of the physical processes 
driving cluster formation and the evolution of cluster galaxies. 

Subsequent studies using CCD detectors have generally 
lacked either the statistics or the radial coverage required 
to test this interpretation. Moreover, most recent measure- 
ments of individual BCG profiles out to large radii (r > 
200 kpc) have f ailed to detect a second component. With 
the exception of iFeldmeier et al.1 (12002ft . these studies find 
profiles that are roughly consisten t with a pure r^/ ^ law 
(Uson.Bou2h. &Kuhn 1990, 1991: Scheic k & Kuhni n^94: 
Gonzalez et al. 2000) . With this program we aim to utilize 
a large sample of BCGs to constrain the ubiquity of the outer 
component and to quantify the properties of this component 
if it exists. 

Our work extends beyond previous studies in several ways. 
First, by utilizing highly uniform drift-scan data we provide 
the first observations that recover the profiles for a statisti- 
cal sample of BCGs to large radii using modern CCD pho- 
tometry. Specifically, we present the surface brightness pro- 
files to r > 300 kpc for 24 BCGs. Second, scanning provides 
long, contiguous images obtained with the same CCD. The 
importance of a precise and accurate sky level is paramount in 
this type of study, and the geometry of drift-scans facilitates 
this endeavor Third, while previous studies typically either 
azimuthally average or model the surface brightness profile 
along the semi-major axis, we employ full two-dimensional 
modelling. This approach enables us to model ellipticity or 
position angle gradients that may be present and to simulta- 
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neously model multiple objects. The latter minimizes the im- 
pact on the profile of other galaxies projected near the center 
oftheBCG. 

The organization of this paper is as follows. In ^we dis- 
cuss the data and preliminary reductions. In ^we describe 
the modelling of the galaxies and present the resulting pro- 
files. In 13 we find that two-component fits provide better 
fits to the profiles than one-component models, compare the 
relative properties (luminosity, effective radius, ellipticity, po- 
sition angle, centroid) of the two components, discuss the fun- 
damental plane for both components, and consider how one- 
component fits impact determination of the BCG total mag- 
nitudes and fundamental plane. We assess the robustness of 
the derived structural parameters in the Appendix. Finally, we 
summarize our conclusions in ^ Throughout the paper we 
assume a cosmology with Hq = 70 km s"' Mpc"', Vl = 0.3, and 
A = 0.7. 

Our results provide important tests for the new genera- 
tion of cosmological cluster simulations that have both the 
force resolution and complexity to generate intracluster stars 
(jSommer-Larsen, Romeo, & Portin ari 2004; M urante et al.j 
|2()04; Willman et al. 2004|) and that are beginning to make 
predictions for the properties of this population. 

2. DATA AND REDUCTIONS 

2.1. Sample Definition 

We observed a total of 30 clusters spanni ng a range of ve- 
locity dispersions and Bautz-Morgan types ( Bautz & Morgan) 
[1270). Our sample consists of nearby clusters (published red- 
shifts 0.08 < z < 0. 12) that contain BCGs with position angles 
within 45 degrees of the east-west axis (our drift scan direc- 
tion). In this paper we focus on the 24 clusters with a single, 
dominant BCG that show no evidence of an ongoing merger 
We note that Bautz-Morgan type III clusters are underrepre- 
sented in the full sample and that errors in the published red- 
shifts yield a wider actual redshift range (0.03 <z < 0. 13). 

2.2. Data 

We obtained drift scan data during four separate runs on the 
Las Campanas Im Swope telescope between 1995 and 2000 
(TableQ. Conditions were photometric during all runs, with 
the data from 2000 having the darkest sky and best seeing. All 
ima ging was obtained in Gunn ; using th e Great Circle Cam- 
era (iZaritskv. Schectman. & Bredthauerill996 ). and the pho- 
tometry is caUbrated to Cousins / using Landolt standards. In- 
dividual drift scans are 2048 x 10000 pixels with a plate scale 
of 0.7" pixel"', except for the scans of Abell 1651, which are 
13000 pixels long. Exposure times per scan vary depending 
upon the declination of the target, ranging from 95-133 sec- 
onds, with seven scans obtained for each target when possible. 
Observational details for each field, including total exposure 
time and seeing, are listed in Tabled For the 24 clusters that 
are the focus of this paper we also include the Two Micron 
All Sky Survey (2MASS) Extended Source Catalog designa- 
tion for the brightest cluster galaxy. 

2.3. Reductions 

The two factors that typically limit photometric measure- 
ments of low surface brightness features are variations in the 
detector response (i.e. flatness) and variations in the sky flux. 
Minimizing residual variations from both sources is criticial 
for this program, as brightest cluster galaxies are expected 
to have /i/ k, 26 mag arcsec"^ at r = 300 kpc. To reliably 



model the profiles out to this radius, we require that this sur- 
face brightness exceed the 3(T systematic uncertainty in the 
sky level. 

2.3.1. Flatfielding 

The basic reduction procedure is similar in form to that de- 
scribed by Gonzalez et al. (2000). With drift-scan imaging, 
pixel-to-pixel variation is mitigated because data are clocked 
across the chip, and so sensitivity variations are a concern 
only perpendicular to the readout direction (at a level of ^ 2% 
in our raw data). After subtracting the bias, we use the scans 
to generate a one-dimensional flatfield, using sigma-clipped 
averaging to compute the sensitivity for each column. The ac- 
curacy of this initial flatfield is limited by both contamination 
from astronomical objects and temporal sky variations along 
the direction of readout. 

We generate a second-pass flatness correction by masking 
objects and modelling the sky variations (details are given in 
the discussion of sky subtraction below). This improved flat- 
field is then applied to the original bias-subtracted images. We 
generate a single flatfield for each of the 1995 and 1999 ob- 
serving runs, while for the 2000 run we have sufficient data 
to divide the run into three two-night subsets and generate 
an individual flatfield for each subset. The difference among 
these flatfields reflects the remaining uncorrected sensitivity 
variations. The rms difference between the second and third 
subsets of the 2000 data is 0.05%, with no detectable sys- 
tematic variation. The flatfield from the first two nights ex- 
hibits a roughly linear systematic variation across the detector 
relative to the other subsets at the level of ±0.17%, with an 
rms difference of 0.1 1%. We conclude that the residual flat- 
ness variations for a single scan in our data are < 0.20% (or 
/i/ ~ 26.4 mag arcsec"^), with column-to-column rms varia- 
tions ~ 0.05% (iij ~ 27.9 mag arcsec"^). Residual variations 
are further mitigated by coadding the scans, and any remain- 
ing gradients (such as the linear difference between flatfields 
noted above) are muted further by the two-dimensional sky 
subtraction procedure described below. 

2.3.2. Subtraction of Sky Variations: First Pass 

The mean sky level typically changes by ^ 3% during sin- 
gle scan. To correct for this temporal variability, we first 
mask all objects in a scan to prevent them from biasing our 
model of the variation. We create a binary mask using seg- 
mentation images output by SExtractor2.3 (Bertin & Arnout^ 
Il996h . excluding all pixels within the detection and buffer re- 
gions around each object. For objects with an area in excess 
of 250 sq. arcsec (primarily saturated stars), the buffer region 
consists of pixels within 150" of the object detection region. 
For smaller objects the buffer radius is either 50" (for those 
objects with area >25 sq. arcsec) or 5" (for those with area 
<25 sq. arcsec). As an added precaution, the entire region 
within 350" of the BCG is masked. 

We then model sky variation along the scan with a five seg- 
ment cubic spline. Subtraction of this model provides a first- 
pass removal of the large-scale temporal variations, which is 
a necessary prerequisite for co-adding. A more precise, two- 
dimensional modeling of the sky is performed at a later stage. 

2.3.3. PSF Modelling 

Following bias subtraction, flatfielding, and the first-pass 
sky subtraction, we register and average the scans to generate 
a single image for each cluster. A key concern at this stage 
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Fig. 1 . — Top: Empirical model PSFs of saturated stars generated for each 
observing run, normalized to the same suiface brightness at 50". The PSF 
for the Mai'ch 1995 ran is the least well determined, being constructed only 
from stars in the field of Abell 1651. We only plot r > 10", where the profiles 
are independent of the mean seeing during the ran. Bottom: Differences in 
the PSFs relative to the November 1999 run. The PSFs for the different runs 
agree to within roughly 30% out to 200" . 



is that the point spread function (PSF) of a bright, saturated 
star can introduce appreciable gradients in the sky level even 
at distances of several arcminutes. To minimize the impact 
of these stars on our photometry, we construct an azimuthally 
averaged model of the PSF out to r = 1000 pixels (700"). A 
single model is generated for each observing run via a lumi- 
nosity weighted combination of all stars with mj < 12.2 in 
images from that run, masking nearby objects and the satu- 
rated cores. The resulting PSF models are shown in Figure 
[2 normalized to have the same surface brightness at 50" ra- 
dius. The PSF for the 1995 observing run is the least well 
constrained because it is generated using only stars from the 
field of Abell 1651. Nonetheless, all profiles agree to within 
^ 30% out to 200", and the difference between profiles at 
700" is not a dominant source of uncertainty in recovering 
the profile of the BCG unless a star with m/ < 7 lies within 
this distance. Abell 2376 and Abell 3809 are the only clus- 
ters in our sample with such a luminous nearby star, while for 
Abell 2405 a bright star lies just outside this radius. 

2.3.4. Subtraction of Sky Variations: Second Pass 

After accounting for the PSFs of bright stars, we proceed 
with a second pass to model the sky. The remaining large- 
scale gradients in the sky level may arise from both instru- 
mental (uneven illumination) and astrophysical (dust emis- 
sion) sources. We construct a background map via wavelet 
decomposition of the input image using the algorithm wvde- 
comp ( Vikhlinin et al. 1998). We use an a trous wavelet trans- 
form to detect all structures with > 2g s ignificance on scales 
of 2° -2'' of 0.7" -358"; see lStarck & Murtagh. 1994 and ref- 
erences therein for a detailed description of the a trous trans- 
form). We subtract these structures from the input image im- 
age and then generate the background map by convolving the 
residual image with a kernel of scale 2'" pixels (716"). This 
scale is larger than the region within which the BCG pro- 
files are modelled (typically r < 250"). We find no evidence 



of oversubtraction in the region of the BCG at the level of 
^, = 27.3 mag arcsec"^ in stacked background images. The 
rms amplitude of the variations removed by subtraction of the 
background maps is equivalent to /^ ~ 26.8 mag arcsec"^. 

Figure |2l illustrates the impact of PSF subtraction and 
wavelet removal of large scale gradients. The three masked 
images are shown at the same contrast and are Gaussian 
smoothed to enhance low surface brightness features. 

3. MODELLING THE SURFACE BRIGHTNESS PROFILES 

We model the brightest cluster galaxy surface brightness 
profiles using GALFIT, a two-dimensional galaxy-fitting pro- 
gram (Peng et al. 2002) that has several advantages for our 
work. First, it permits simultaneous modelling of multiple 
components. This feature enables us to model galaxies near 
the center of the BCG that could otherwise bias the derived 
profiles, to model systems with multiple dominant galaxies 
(in a subsequent study), and to model the BCGs themselves 
with multiple components. Second, because we fit the two- 
dimensional luminosity distribution rather than the surface 
brightness profile along the major axis, the x^ of the fit is 
directly sensitive to variations in the position angle and el- 
lipticity with radius — unlike many previous studies. The 
gradients provide important leverage in discerning wh ether 
two compo nent models for BCGs (e.g ., ISchom bert 198^ 
fl988: .Porter. Schneider. & Hoesseil991h are superior to sin- 
gle component Sersic and r'/'* models. Finally, the program 
generates an image of the best-fit model, from which we cal- 
culate the model's radial profile using the same masking as 
used for the data. This radial profile is valuable both for visu- 
alization and comparison with previous studies. 

There are five critical inputs to GALFIT, aside from the im- 
age itself, that warrant detailed discussion. These inputs are 
the stellar PSF, the image mask, the noise map, the mean sky 
level, and the parametric model for the surface brightness dis- 
tribution. Each is discussed below, and the sensitivity of the 
results to the PSF, masking, and sky level are explored in the 
Appendix. 

3.1. PSF 

GALFIT requires a representation of the stellar PSF for its 
convolution with the galaxy model. We der ive the PSF us ing 
the IRAF implementation of DAOPHOT II JStetsonll987h . A 
Moffat model with first order residual variations is generated 
for each image using the 150 brightest unsaturated stars. 

3.2. Masks 

Masks are generated using the procedure described in |2] 
A key difference is that here we do not automatically mask 
any objects within 20" of the center of the BCG. We visu- 
ally inspect the images, manually masking these objects when 
possible and leaving them unmasked for modelling when nec- 
essary. We adopt the latter option only when part of the ob- 
ject Hes within - 5" of the BCG. Two clusters, Abell 2969 
and Abell 3705 require special attention. In Abell 2969 two 
galaxies overlapping the BCG were modelled, while in Abell 
3705 a single large, bright galaxy 2' away was modelled to 
minimize its impact upon the BCG profile at large radii. 

For each BCG, we also exclude all pixels outside a circular 
aperture of physical radius 300 kpc. We choose a fixed phys- 
ical radius to facilitate comparison of the profiles, and choose 
300 kpc because this is roughly the distance at which the sur- 
face brightness along the semi-major axis is comparable to the 
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Fig. 2. — Series of images of Abell 0122 illustrating the impact of PSF subtraction and wavelet removal on the level of sky variation. The size of the region 
shown is 30' X 19'. The top left image is the data immediately after co-adding, with the saturation level corresponding to /.j = 21.8 mag arcsec"^. Subsequent 
images are Gaussian smoothed to enhance the visibility of low surface brightness features and are each displayed with the same contrast. The saturation level in 
these images corresponds to /x = 25.8 mag arcsec"^ . We mask objects other than the BCG, which is marked by an overlaid circle of radius 300 kpc (the size of 
the fitting region in ^. The upper right image, like the upper left panel, is constructed from the coadded individual scans. Significant excess flux is visible along 
the perimeters of stellar masks. The lower left panel demonstrates the improvement achieved with the PSF subtraction of saturated stars. Finally, the lower right 
panel demonstrates the improvement achieved, particularly in the uniformity of the background on large scales, by subtracting the wavelet background map. 



5(7 sky level uncertainty for most BCGs in our sample. For a 
typical BCG, ^25-45% of pixels within this physical aper- 
ture are masked. In only three cases does the masking exceed 
50% of the pixels within this radius: Abell 2376, Abell 2405, 
and Abell 3705, with ^ 55% of pixels masked for each. For 
Abell 2376 and Abell 2405 the increased masking is due to a 
nearby bright star; for Abell 3705 the cluster simply lies in a 
region of high stellar density. 

3.3. Noise maps 

Noise maps are constructed directly from the data, which 
has zero mean sky due to previous processing. Objects are 
initially masked and the DIMSUM^ routine iterstat is used to 
directly compute the background rms noise level. Assuming 
Poisson statistics, we add the equivalent sky level to the im- 
age and take the square root to generate a noise map, which 
is then smoothed with a 2 pixel FWHM Gaussian. Smoothing 
provides a more robust measure of the noise level, and use of 
a smoothing scale comparable to the seeing minimizes spatial 
degradation of the noise map. While these noise maps yield 
model fits that are good in absolute terms (reduced x^ ~ 1; 
©, because these noise maps are approximate we only draw 
conclusions from the relative x^ values derived for the differ- 
ent models. 

3.4. Sky Level 

The background level is a critical input because systematic 
error in the sky will dominate the profile at large radii. We use 

' DIMSUM is the Deep Infrared Mosaicing Software package developed 
by Peter Eisenhardt, Mark Dickinson, Adam Stanford, and John Ward. 



GALFIT to determine the sky level, utilizing the same noise 
maps and masks as for the profile fitting. The sky level is mea- 
sured within circular annuli of radii r = 600-750 kpc centered 
on the BCG. The extrapolated flux level of the central galaxy 
at these radii is more than 3 mag fainter than it is at 300 kpc 
(or equivalently <7% as bright), which is 2-3 times less than 
our Icr sky level uncertainty (see below), and hence negligibly 
impacts our estimate of the sky level. 

We use the ensemble properties of our sample to quantify 
the uncertainty in the derived mean sky levels. The measured 
dispersion among the mean sky levels observed for different 
galaxies {aohs) has two sources: true differences in the global 
mean sky levels of different images {(jg) and our measurement 
uncertainty within the background annuli (a„,). If we assume 
that these contributions add in quadrature, then ct,^, = T^ij^-cr^. 
To estimate cjg we measure the global sky level over the entire 
image (excluding the BCG region) for each cluster and then 
compute the dispersion. We then measure the mean sky level 
within the background annulus for each BCG and compute 
Gobs- In both cases we use GALFIT to determine the mean sky 
level, and together these data enable us to derive (t,„, the \a 
uncertainty in the sky level. For a typical galaxy in our sample 
this 1(7 uncertainty corresponds to /i ci 28.4 mag arcsec"^. '^ 
In one case, Abell 2376, we increase our estimate of the sky 
uncertainty by a factor of two to account for the impact of a 
proximate bright star. 

3.5. Input Models 

■* The precise value varies depending on the photometric zeropoint of the 
field, because the calculation is performed using uncahbrated counts. 
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We test three different functional forms for the BCG surface 
brightness distribution. Two of these are single component 
deVaucouleurs (r'/^) and Sersic (r'/") models, while the third, 
in the spirit of Schombert ( 1988), is a two component model 
consisting of a pair of r'/^ profiles with independent param- 
eters. We will refer to the latter model as the 2-deV model 
henceforth. The r'/" model has one more degree of freedom 
(dof) than the r'/'* model, while the second component in the 
2-deV model adds six dof relative to the single r^l^ model 
(x,y,^e,re,0,e). For each functional model we run GALFIT 
with a series of different initial parameters to ensure conver- 
gence to the global minimum for each model, and iterate upon 
the solution until Ax^ < lO""*. 



4.1.1. Sersic versus r '' '^ 

Our data demonstrate that the r'/" profile provides a dra- 
matic improvement relative to the r'/'' profile. The addi- 
tional degree of freedom decreases x^ by an average of 3649 
(Table O, and substantially improves the fit at large radii 
(Figures |3]— E)- In only two cases is the decrease in Ax^ 
minimal (Abell 1651 with Ax^ = 38 and Abell 2721 with 
Ax^ = 19-84 for data from two different years). Abell 1651 
is the cluster from our original pilot study in which we found 
that the BC G showed no strong evidence of departure from 
r'/'* (Gonzalez etal."2000). We Hst in Table|5lthe Sersic pa- 
rameters derived in both studies of this BCG, noting that the 
results are consistent despite different modelling procedures 
and bandpasses. 



4. RESULTS AND DISCUSSION 
4.1. Evidence for Two Components 

The optimal parameters for each model and the Ax^ val- 
ues between different models are presented in Tables |2l— 
In these Tables we include both the statistical parameter un- 
certainties from GALFIT and, in parentheses, the uncertain- 
ties corresponding to Icr changes in the sky level (see ^3.41 
and the Appendix). To illustrate the quality of the models, 
we also present one-dimensional profiles extracted using the 
IRAF task ellipse. In Figure|3]we present surface brightness 
profiles along the semi-major axis, and in Figure |3 we show 
the residuals relative to the 2-deV model. In both figures we 
employ the same masking and elliptical apertures for the data 
and models to enable direct comparison.'' The model images 
output by GALFIT include convolution with the input PSF, 
and hence comparison to the data is valid over the entire radial 
range shown. In Figure |4] we overlay shaded regions corre- 
sponding to the 1(7 systematic sky level uncertainty (0^^3) 
to clarify the significance of deviations between the data and 
models. 

Figures |5] and |6l show the ellipticities and position angles 
derived for both the data and best-fit models. While the sin- 
gle component models have fixed values for both quantities, 
there are gradients visible in the Figures for these models be- 
cause we have folded into the profiles several additional fac- 
tors to enable direct comparison with the data. First, masking 
of galaxies near the center of the BCG can lead to variations 
in both quantities. We use the same masking for both the data 
and models to ensure a fair comparison. Second, seeing acts 
to decrease the central ellipticity and can cause position angle 
changes if the ellipticity is small. The models are convolved 
with the PSF from the image and so exhibit the same decrease 
as the data. Third, in some instances we found it necessary 
to model rather than mask secondary galaxies near the center 
of the BCG, which are listed in the caption for Figure |5] The 
flux from these galaxies is included in the plotted profiles for 
both the data and models. Abell 2400 is a good example of the 
impact that such a companion can have on the central profile. 
Together the above factors explain the gradients seen in the 
single component models; yet, even including these factors, 
these models cannot reproduce the ellipticity and position an- 
gle variations observed in some of the BCGs. 

^ A potential concern is that the lack of ellipticity and PA gradients for the 
Sersic and r' " models could bias the residuals shown in Figure|4] We have 
performed tests to assess the importance of this bias, and find that the effect 
is negligible. 



4.1.2. 2-deV versus Sersic 

Although the Sersic profile is superior to the r'/"^ model for 
our BCG sample and provides a fit that is good in absolute 
terms (xj, — 1), the Ax values in Tablel^argue that the 2-deV 
profile provides the best parametric fit to the data. With five 
additional degrees of freedom, the 2-de V profile decreases x^ 
by an average 1039 relative to the Sersic profile, which is 
highly significant. The reasons that the 2-deV model provides 
a superior fit can be assessed by inspection of Figures |3l— El 
For roughly a third of the galaxies the two components are 
required to reproduce the observed large gradients in the el- 
lipticity and/or position angle. In the other two-thirds, where 
these gradients are small, the fit is still statistically superior. 
The 2-deV model is thus preferred due to both its ability to 
fit the mean surface brightness profile and improved ability to 
reproduce the gradients. 

4.1.3. Physical interpretation of 2-deV model 

Does the success of the 2-deV model imply that there are 
two physically distinct components? Perhaps the ellipticity 
gradients and isophotal twists are simply the result of viewing 
triaxial galaxies off-axis. If so, then the individual compo- 
nents are not physical but rather a convenient means of ap- 
proximating a triaxial stellar distribution. There are several 
reasons to conclude that this is not the case. First, there are 
radial regions over which the position angle and ellipticity are 
constant (e.g. Abell 2571), which implies that at those radii 
we are either viewing the system along a principal axis or else 
that the system is not triaxial (see Figures^— |6|l. If the former, 
then the ellipticity and position angle should remain constant 
at all radii. Second, the more dramatic position angle shifts 
(e.g. Abell 2571, Abell 3166) happen very suddenly, which 
again should not be the case for a triaxial system. Third, the 
apparent position angle of a triaxial body is in general differ- 
ent than the projected position angle of the true major axis. 
The kno wn corr elation between BCG and cluster orientation 
( Bing^lj[l982) argues against these galaxies being strongly 
triaxial systems viewed at random angles. Fourth, if we are 
viewing projected triaxial ellipsoids, then we should expect 
all the galaxies to show position angle swings, as is true with 
most ellipticals. Many of the BCGs in this sample in fact 
lack significant position angle gradients, and it is statistically 
quite improbable that we are viewing such a large fraction 
along their principal axes. Finally, for a triaxial system the 
relative importance of the apparent second component should 
be a strong function of viewing angle (i.e. when the BCG 
is observed near the principal axis this component should be 
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Fig. 3. — /—band surface brightness profiles for the 24BCGs modelled in this paper. The error bars, which are typically smaller than the symbols, are statistical. 
Overlaid are the corresponding surface brightness profiles for the best-fit two-dimensional models. We emphasize that these are not fits to the one-dimensional 
profile, but rather measurements of the two-dimensional models taken in the same apertures as for the data. The r'''' profile (short-long dashed) is a poor 
representation of the data in most cases, while the r'^" (dotted) and 2-deV (solid) profiles provide improved fits. The solid and dashed horizontal lines correspond 
to the 3(T and 5o" sky level uncertainties, respectively. 
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Fig. 3. — Continued. 
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Fig. 4. — Surface brightness residuals relative to the 2-deV model. The line types are the same as in Figure|3] but we have now also added a shaded region 
denoting the Icr systematic sky level uncertainty. We note that in two cases, Abell 2721 and Abell 3705, bright saturated stars appear to skew the residuals at 
large radii. In both cases the stellar PSFs were subtracted and the central regions masked, but the systematic sky level uncertainty is likely larger than our estimate 
in these cases. 



negligible). We find instead that the scal e len gths of the two 
components are strongly correlated (see ^4.3t . 

Aside from the above arguments against the observed gra- 
dients in ellipticity and position angle being due to triaxiality. 



ancillary data also favor the two components having physi- 
cal significance. Most compelling are the velocity dispersions 
in two nearby clusters, Abell 2 029 (Dressier 1979) and Abell 
2199 tCarter & Bridgeslll999tlkelson et alJI2002l) . which are 
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Fig. 5. — Ellipdcity profiles for the 24BCGs modelled in this paper. The line types are the same as in Figurey] Note that the sharp, non-monotonic changes in 
ellipticity and position angle seen in a few of the profiles (e.g. Abell 0122) occur when a significant fraction of the data at that radius are masked. These features 
are thus more a function of the masking than indicative of the profiles, but are produced in both the model and data profiles since the same masking is appHed in 
both cases. In contrast, large monotonic changes in the ellipticity (e.g. Abell 3112) are physical in origin. The 2-deV profile is best able to reproduce these large 
monotonic changes. 
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observed to rise towards the cluster velocity dispersion at 
large radii. Such a rise is strong evidence that the orbits of 
the outer stars are associated with the gravitational potential 
of the cluster rather than the central galaxy. Taken with the 
arguments against triaxiality, these data argue that the success 
of the 2-deV model is not simply a result of analyzing triax- 
ial bodies, but rather that these two components are indeed 
distinct physical entities. 

Finally, we emphasize three additional points. First, while 
we model the two components as a pair of r'/'* profiles, we are 
not specifically concluding that the two components are r'/"^ 
in nature. This model simply provides an adequate description 
of our data with a minimum number of free parameters. The 
outer componen t could conceivably be b etter described by an 
NFW profile cNavarro. Frenk. & Whitelll996. .IQQ?.) . but we 
have not yet explored this model. Second, for any single 
BCG, we caution the reader to not place undue significance 
on the structural parameters of the individual components in 
the 2-deV model. While the results are robust for the statis- 
tical sample, the decomposition for a given BCG may be bi- 
ased by internal surface brightness structures or tidal features. 
Third, we are not claiming to see the dramatic profile breaks 
characteristic of classical cD envelopes (e.g. Schombert 198q, 
0588) as ubiquitous features. The transition between the two 
components in our data is typically much more subtle. 

4. 1 .4. Evidence for additional components 

Given that we find a two component model to be statisti- 
cally superior, one can ask whether there is evidence for any 
additional components. The residuals in Figure|4]indicate that 
in a few cases the flux exceeds the best-fit model by several 
sigma in the outermost few radial bins. In two of these cases, 
Abell 2721 and Abell 3705, the likely origin of this excess is 
residual flux from saturated stars at the corresponding radii. 
In another, Abell S0296, only the 2000 observations show an 
excess, indicating that the origin is likely instrumental rather 
than physical. The other two cases with ^ 2ct excesses (Abell 
165 1 and Abell S0084) lack a similar explanation, and we also 
note that the sample as a whole exhibits a slight tendency for 



excess flux at large radii. With the current data, however, we 
cannot discern whether these excesses are due to a systematic 
error or a true physical excess relative to the model, arising 
from either an additional component or the outer profile hav- 
ing a form other than r'/"^ . 

4.2. Relative Properties of the Two Components 

4.2.1. Luminosities 

What fraction of the light does the outer component con- 
tribute? Figure0demonstrates that the outer component typi- 
cally contains ~ 80-90% of the total light from both compo- 
nents. In only five cases is the outer component's luminosity 
<70% of the total, and those systems tend to be the ones with 
the largest uncertainties in the magnitude of the outer compo- 
nent (e.g. Abell 1651, Abell 2400, Abell 2405, APMC 020). 
In subsequent plots we mark these five systems with starred 
symbols to illustrate the impact of this systematic difference 
in luminosity ratio (whether real or due to an error in decou- 
pling the two components). 

4.2.2. Effective Radii 

How do the effective radii of the two components compare 
and are they correlated? The range of effective radii for the 
inner component is 1-48 kpc, while the range for the outer 
component is 18-480 kpc. Figure |8l shows that the effective 
radii of the two components are coupled, with a typical ratio 
fe.oiiier/fejmier ^ 10-40. For Comparison, we overlay a line 
denoting the ratio of the projected half-mass radii for the dark 
matter and stellar components (rDM/fswiiar = 13.2) from the 
cluster simulation of Dubinski ( 1998). This simulation recov- 
ered a purely r'/'* BCG profile, with a scale comparable to 
that of the inner components observed in this study, but the 
stellar component was not modelled with sufficient resolution 
to produce excess light at large radii. 

4.2.3. Fundamental Planes 

Do the individual components in the 2-deV model indepen- 
dently obey a < p > -r,, relation similar to the projection of 
the elliptical galaxy fundamental plane? We find that both 
components do indeed form tight < p> -r^ relations, as seen 
in Figure |9j?- The slope of the relation for the inner compo- 
nent is consistent with the elliptical fundamental plane, which 
has < / >cx r^'^^ ( J0rgensen, Franx, & Kiaergaard 1996). In 
constrast, the slope of the relation for the outer component is 
steeper than for elliptical galaxies. If we exclude systems in 
which the outer componen t contr ibutes less than two-thirds of 
the total luminosity (see ^4.2.11 above), the slope of the re- 
lation is < / >(x rj'-^ and the scatter is consistent with the 
observational uncertainties (Figurel^i)- 

4.2.4. Ellipticities and Position Angles 

Two other parameters of basic interest are the elliptici- 
ties and relative offsets in position angle of the two com- 
ponents. The outer component is generally more elliptical 
than the inner component (Figure^p), consistent with previ- 
ous ob servations that found increasing ellipticity with radius 
(c.f. P orter. Schneider. & Hoessell 19911) . The ellipticity of the 
outer component is comparable to typical observed elhptici- 
ties for the distribution of cluster galaxies ("Carter & Metcalfd 
1980; Plionis. Barrow. & Frenk 1991; Buote & Canizarel 
1996; Basilakos. Phonis. & Maddox 2000). 

The inner and outer components are closely aligned ( A0 < 
10°) for roughly 40% of our sample, indicating that both may 
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retain information about the initial infall direction of the sys- 
tem (FigureFTOb ). Those inner components that are not well- 
aligned show a preference for large A6 (A9 > 45° with a dis- 
tribution that rises toward 90°). This offset may be the result 
of infall events along other directions that determine the cur- 



rent position angle of the central galaxy. It is also possible that 
this could be a simple projection effect in some cases — par- 
ticularly cases in which A9 ^^ 90° and b/a > 0.9 (e.g. Abell 
3727). 
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Fig. 1 1. — Difference in absolute magnitudes between the Sersic and 2- 
deV models. The magnitudes are comparable for 75% of the galaxies, while 
for the other 25% the Sersic magnitudes are >0.5 mag brighter. 



4.2.5. Centroid Offsets 

Finally, are the two components concentric? When mod- 
elling these galaxies with GALFIT we allowed the centroids 
of the inner and outer components to vary independently. We 
find that the two components are concentric to within 1.2" 
(less than the seeing FWHM) in all but one case. The ex- 
ception is Abell 1651, in which the centroids are offset by 9" 
(- 15 kpc). 

4.3. Biases Arising from One-Component Models 

In preceding sections we discerned that the 2-deV model 
provides a better representation of the typical BCG surface 
brightness distribution than do r'/"* or r'/" profiles, and quan- 
tified the properties of the components in this model. We now 




100 



Fig. 12. — Surface brightness vs. effective radius for the single r'" pa- 
rameterization. The slope of the solid line corresponds to the < fj, > —re 
relation from J0rgensen, Franx. & Kiaereaard 1 1996), assuming constant ve- 
locity dispersion and arbitrarily normalized to match the present data. Despite 
the failure of the r'''' model at large radii, this parameterization yields a tight 
< fi > —Ke relation with a slope similar to the relation for normal elUptical 
galaxies. 



return to the single component models, both to facilitate com- 
parison with previous work and to explore the relative biases 
that arise from using these different models. 

4.3. 1 . Effects on Total Magnitude 

How does the choice of parametric model influence the total 
magnitude derived for the brightest cluster galaxy? The bulk 
of published BCG photometry assumes either an r'/"^ or r'/" 
profile to calculate a total magnitude. The r'/"^ model obvi- 
ously underestimates the total luminosity relative to the Sersic 
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Fig. 13. — Parameters for the Sersic model, (a) Sersic index vs. effective radius. The dashed arrow illustrates the intrinsic correlation between re and 1/n in 
the Sersic model. To generate this correlation vector we fit models to a single galaxy multiple times with n fixed to different values. The data indicate that the 
true uncertainties in the Sersic parameters are dominated by degeneracies among the con'elated parameters and are larger than the quoted uncertainties. Note that 
the length of the correlation vector is not meant to correspond to the size of the con 'elat ed uncertainties, but is simply set to span the range of observed re values. 
(b) Surface brightness vs. effective radius. The solid line is the same as in Figure fT2l while the dashed arrow again indicates the correlation vector associated 
with changes in n. For large effective radii the data lie parallel to this correlation vector. It is thus difficult to decouple any hue deviation from the FP relation 
from correlated errors. 



and 2-deV models, while Figure |3] argues that in most cases 
the Sersic and 2-deV models should yield similar total mag- 
nitudes. Figure ^2 demonstrates that this expectation is cor- 
rect for roughly three quarters of the galaxies. For the other 
25%, the magnitudes derived from the Sersic parameteriza- 
tion are more than half a magnitude brighter (albeit typically 
with large uncertainties). In these cases, the optimal Sersic 
models have n>9 and thus contain significant luminosity at 
r > 300 kpc (i.e. outside the fitting region). 

Table |6l lists the ensemble average and dispersion in the 
total magnitudes derived using each parameterization. We 
also include in this table the magnitudes within circular met- 
ric apertures of radii 300, 50, and 10 kpc. Within the metric 
apertures the Sersic and 2-deV models yield consistent mag- 
nitudes and dispersions, while the r'/"* model underestimates 
the flux by ^ 30% within the 300 kpc aperture. For the total 
magnitudes, the Sersic model yields significantly larger scat- 
ter, which is a consequence of the extrapolation at r > 300 kpc 
for the galaxies with large n described above. We note that 
the dispersions for the smallest apertures are consistent with 
[ Postman & Lauer ( 1995), who found a dispersion of 0.327 
mag within a metric aperture of 10 /jJqq kpc. 

4.3.2. Effects on the Fundamental Plane 

What is the effect of fitting a one-component model on 
the fundamental plane derived for BCGs? Figure [T2l shows 
the < fi > -Te relation for the single de Vaucouleurs model, 
which has a slope similar to that of the < /i > -r^ projec- 
tion of the fundamental plane (FP) for normal elliptical galax- 
ies (|j0rgensen. Franx, & Kiaergaard 1996). It is interesting to 
note that the galaxies roughly obey a fundamental plane type 
relation for the r'/"* fit despite the demonstrated inferiority of 
this model. Figure[T^shows the l/n-r^ and < /i > -r^ rela- 
tions for the Sersic model. Because all three parameters are 
strongly correlated, we explore the impact of changing n on 



the other parameters. The additional degree of freedom pro- 
vided by n yields much larger effective radii (the largest r^ 
in the r'/"* fits is < 60 kpc) and fainter surface brightnesses 
(the faintest < fi > in the r'/"* fits is < 23). The slope of the 
< fi> -Ke relation is comparable to the direction of the over- 
laid correlation vector at large r^ (i.e. large n), making any 
true deviation from the FP relation difficult to decouple from 
correlated errors arising from uncertainty in n. The scatter 
in the 1 /n - r^, plot is also predominantly along the direction 
of the correlation vector Together these plots illustrate a key 
limitation of the Sersic model for this analysis — namely that 
systematic variations in n can induce large correlated uncer- 
tainties in the other parameters. 

5. CONCLUSIONS 

Using highly uniform photometric data, we study the sur- 
face brightness distribution of the brightest cluster galaxy 
(BCG) in each of 30 clusters at 0.03 < z < 0.13 that span a 
range of velocity dispersions and Bautz-Morgan types. In this 
first paper, we focus on the 24 clusters for which there is a sin- 
gle, dominant BCG and no evidence of a substantial ongoing 
merger of subclusters. We employ full two-dimensional fit- 
ting to model the surface brightness profiles of these galaxies 
and their stellar halos within r = 300 kpc. From this analy- 
sis, we reach several conclusions regarding the nature of the 
diffuse optical light in clusters: 

1. The observed profiles are well described by a fwo- 
componenf model consisting of two r'/^ profiles with different 
scale lengths, ellipticities, and orientations. This model yields 
significantly lower x^ than single r^l* or r'/" profiles and re- 
produces the position angle and ellipticity gradients observed 
in these galaxies. The nature of these gradients is such that 
they are unlikely to originate fro m v iewing a single compo- 
nent triaxial system at random ( ^4.1> . We conclude that the 
two components should be considered physically distinct. 
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2. The inner component is clearly associated with the BCG. 
This component is structurally similar to a normal bright el- 
liptical galaxy. The slope of the < /i > -r^ projection of the 
fundamental plane is consistent with published results over 
a comparable range of radii (r^. « 1 - 30 kpc), and the ob- 
served axis ratios {b/a = 0.77 ± 0. 14) are comparable to those 
obser ved in normal ellipticals (iPorter. Schneider. & Hoessell 
Il99l . 

3. The outer component is consistent with arising from a 
population of intracluster stars. This intracluster light (ICL) 
component dominates the luminosity, containing roughly 80- 
90% of the total light of the two components. The large physi- 
cal scale of the ICL component (r^ 2± 20-500kpc) is similar to 
that expected for the dark matter distribution (Dubinski 1998) 
and is thus consistent with the picture in which the intracluster 
stars trace the global cluster potential. The ICL component 
exhibits a tight < fi > -r^ relation, but with a steeper slope 
than the corresponding projection of the elliptical galaxy fun- 
damental plane (< / >cx r~^^^ compared to < / >oc rj'^', re- 
spectively). The ICL component is typically more elongated 
than either the inner BCG component or normal ellipticals, 
with axis ratios comparab le to thos e observed for the distri- 
bution of cluster galaxies (*BinggelLi .l982) . The existence of 
a distinct ICL component associated with the cluster poten- 
tial is further bolstered by kinematic data that detect rising 
velocity dispersions at r > 10 kpc in several systems (Abell 
2029, Dressier ( 1979); Abell 2199, Carter & Bridges ( 199^; 
iKelson et al.. (,2002)), presumably due to the transition from 
the bound stellar population of the BCG to the intracluster 
stellar population. 

4. The inner and outer components are strongly aligned 
(I A6'| < 10°) in roughly 40% of the clusters. When they are 
not aligned, the components tend toward high degrees of mis- 
alignment, suggesting that accretion of infalling material may 
change the orientation of some BCGs for a time. 

5. The ratio of scale-lengths for the outer and inner com- 
ponents, re^outer/re.inner, IS ~ 10-40, with the distribution 
strongly peaked at r,,outerl runner = 35. 

Recent hydrodynamic simulations incorporate much of the 
relevant physics necessary to produce populations of intr- 
acluster stars, and thus provide quantitative predictions for 
the properties of the ICL. These simulations agree on sev- 
eral interesting points. Consistent with pre vious theoreti - 
cal work (c.f. [Richstone & Malamuth 19831 lMilledll983l) . 
the simulations predict that intracluster stars produce roughly 
20-50% of the total cluster lumin osity and orig inate from 
tidal stripping ( Sommer-Larsen. Romeo. & Portinaril l2004t 
iMurante et a'Lri2004: .Willman et al.. .2004) . The intraclus- 
ter stars in these simulations form at higher median red- 
shift than stars that remain in cluster galaxies at z = 0, 
concentrate more toward the cluster center than the galax- 
ies, and exceed the stell ar mass density w ithin about half 
the virial radius ( ^ Sommer-Larsen. Romeo. & Portinari 2004i 
IMurante et al.l2004l) . Additional! v. formation of the ICLcom- 
ponent is found to be a continuing process, with significant 



growth in luminosity since z=l dWi llman et alJ2004 l). 

We confirm that the ICL is a ubiquitous feature for clusters 
with a single, dominant BCG, demonstrate that it dominates 
the combined luminosity of the BCGh-ICL components, and 
quanti fy some o f the properties of the ICL. Recent work by 
Feldmeier et alJ (i2004 ) indicates that intracluster light is also 
prevalent in clusters without dominant central galaxies, and 
thus together these studies confirm that ICL is a generic fea- 
ture of the cluster environment. Our measurements of the ex- 
tent of the ICL component and its elongation, which is similar 
to published cluster galaxy spatial distributions, further indi- 
cate that the evolution of the ICL is tied to the cluster as a 
whole rather than to the BCG. 

Finally, we emphasize that understanding the growth and 
accretion of brightest cluster galaxies and the intraclus- 
ter stellar populations is also of direct relevance to stud- 
ies of lower mass elliptical galaxies. Graham etal. (199a) 
interpreted the large Sersic indices observed at r < 25 
kpc for their BCG sample as the extension of an em- 
pirical correlation between n and galaxy mass observed 
by Caon. Capacc ioli. & D 'Onofrio ( 1993) for other elliptical 
galaxies, and Graham. Er win. Caon. & Truiillo (2001) have 
presented evidence for a correlation between n and the cen- 
tral black hole mass. We have now demonstrated that simi- 
larly large Sersic indices are obtained when BCGs are mod- 
elled to much larger radii. Because the outer component in 
the 2-deV model contributes a significant fraction (>35%) of 
the total luminosity within r = 25 kpc for 75% of the BCGs 
in our sample, we argue that these large Sersic indices are 
better interpreted as arising from the dominance of this intr- 
acluster component at large radii. It is intriguing to consider 
whether the correlations observed for lower mass systems are 
partly due to a correlation between galaxy mass and the frac- 
tional luminosity existing in halos comprised of stars tidally 
stripped from lower mass satellites. We suspect that investi- 
gation along these lines may provide a physical motivation for 
the observed empirical correlations. 
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APPENDIX 

ROBUSTNESS OF RESULTS 

We now explore whether either our abiUty to discriminate between models or the resulting model parameters are compromised 
by systematic errors and compare results for two clusters observed in multiple runs. 

Error in the sky level 

While GALFIT returns statistical uncertainties, we are undoubtedly limited by systematic errors in deriving structural param- 
eters. The most significant source of such errors is the determined mean sky levels. In this section we ask two questions: 1) 
does systematic error in the sky level degrade our ability to discriminate between Sersic and 2-deV models, and 2) what is the 
magnitude of the induced uncertainties in the structural parameters? 

We estimated the 1 a uncertainty in the background sky levels in ^3.41 To assess whether this degree of uncertainty compromises 
our ability to discriminate between models, we rerun GALFIT for all models with the sky level artificially altered by ±1ct. We 
find that the resulting y^ differences between the r'/" and 2-deV models are Ax^ = 1039 and Ax^ = 1089 when we increase 
and decrease the sky, respectively, which are consistent with our result from ^(Ax^ = 1039). We conclude that the statistical 
preference for 2-deV models is not a result of systematic error in the sky level. 

From the same artifically altered data, we rederive the structural parameters. The errors quoted in parentheses in Tables 12 141 
correspond to the variations arising from this change. For the Sersic model the most sensitive parameters are r^, and «, because a 
change in the outer structure can mimic a change in the sky level. The derived ef fective radii and surface brightnesses are relatively 
robust for cases with n < 10, but are highly uncertain for larger indices (see ^4.3.2> . For the 2-deV model the parameters are 
remarkably robust. The largest uncertainties arise for the outer component (for the same reason as with the Sersic model), and 
the uncertainty is proportional to the effective radius of the outer component. Still, even for the cases where r^ > 300 kpc (i.e. 
the size of the fitting region) the uncertainty in r^ is < 10%. These cases should obviously be viewed with caution (as is always 
true when the fitting region is < r^), but this robustness is encouraging. 

Sensitivity to input PSF 

Another concern is whether the results are sensitive to the accuracy of the input PSF. We assess this sensitivity by reanalyzing 
one of the BCGs (Abell 0122) using PSFs taken from other fields with seeing comparable to within 0.1" (Abell 2376, Abell 
3727, and Abell S1096). Although drawn from similar seeing conditions, these PSFs are structurally distinct. The optimal 2-deV 
parameters are listed in Table0for each input PSF. The most sensitive parameter is r^ for the inner component, which increases 
up to 25% when an incorrect PSF is used. We consider 25% to be an approximate upper bound on the systematic uncertainty in 
the inner r^ due to error in the input PSF. The outer component re, and < fi > for both components, agree to within 7% in all 
cases. Finally, x^ changes by < 100, which is significantly less than the differences between the three functional forms tested in 
this paper. 

A different concern is that the PSF convolved with the GALFIT models only extends to r = 15". Could the inclusion of the PSF 
wings (Figure[3 affect the results? Convolution of the BCG profile with the full PSF including the extended wings increases the 
surface brightness by ^ 1% at r > 10", which should be small enough to not significantly alter the derived structural parameters. 
To verify this assumption, we performed modelling tests employing a large PSF that includes these wings. We find that the 
most sensitive model parameter is the r^ of the outer component in the 2-deV model, which typically decreases by 5- 10% 
when the extended PSF is used. The relative luminosities of the two components also change slightly, with the outer component 
contributing ~ 2% less of the total luminosity. In both cases this bias due to using the smaller PSF corresponds to < Icr changes 
and thus will not qualitatively alter any of the results described below. 

Sensitivity to masking 

Might the "second component" observed at large radii be due to residual emission from individual cluster galaxies? If the 
masking is insufficient, then flux from the outer regions of individual galaxies will appear as a secondary halo in the surface 
brightness modelling. This component would have an ellipticity and orientation similar to the galaxy distribution, potentially 
yielding gradients in both quantities as observed in our data. 

We test this possibility by varying varying the degree of masking and measuring the difference in the observed profiles. We 
perform this test on the BCG in Abell 2571, a galaxy for which the data indicates a significant excess relative to r'/"^ at r > 20 
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kpc. We extend the object masks and increase the fraction of masked pixels from 26% to 53%. The best-fit parameters change by 
less than the systematic uncertainties derived for changes in the sky level, and the 2-deV model maintains a significantly lower 
X^ than the other models ( X^ev ~ ^2-deV =4300 and Xsersw~X2-deV ~ 1745 ). The profiles derived using the normal and enhanced 
masking are also consistent to within the la statistical uncertainties for all radial bins. We conclude that the observed outer 
component is not an artifact arising from insufficient masking of cluster galaxies. 

Clusters observed during multiple runs 

The two clusters in our sample that were observed in both 1999 and 2000, Abell 2721 and Abell S0296, provide an independent 
check on the impact of systematics. We model the BCGs independently for each data set and compare the profiles and derived 
parameters for the BCGs. For both galaxies the azimuthally averaged surface brightness profiles are consistent to within the 
statistical uncertainties at all radii larger than the seeing. With only two exceptions, the parameters derived from the two data 
sets are consistent to within the systematic uncertainties for all models (Table |2]—|4}. The only exceptions are the ellipticities 
for the Sersic and 2-deV models and the position angle shift, A^, in the 2-deV model, which differ in some cases by as much 
~ 5(7. We thus caution that the precision (typically ^ 0.01 in ellipticity and ~ few degrees in position angle) may be somewhat 
overestimated. 

The most striking difference is in the Ax^ values for Abell S0296, which are significantly smaller for the 1999 data set. The 
likely origin of this discrepancy is the poor data quality for the 1999 observations of this cluster, which, as noted in Tabled are 
compromised by large-scale sky variations. Consequently, the comparison for this cluster can be viewed as a worst-case test. 
Nonetheless, both data sets yield the same qualitative conclusions and comparable structural parameters, which again argues that 
our results are robust against systematic errors. 
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TABLE 1 














Cluster Sample 








BM 


Date 


Exposure 






2MASS 


Cluster 


ZBca 


type 


of Run 


time (s) 


Airmass 


FWHM(") 


ID Comments 


Abell 0122 


.1127 


I 


10/2000 


635 


1.03 


1.38 


2MASX J00572288-2616528 


Abell 0447 


.1123 


I 


10/2000 


667 


1.10 


1.67 


a b 


Abell 1651 


.0853 


I-II 


3/1995 


285 


1.12 


1.47 


2MASX J12592251-0411460 


Abell 2376 


.0891 


I-II 


10/2000 


673 


1.19 


1.33 


2MASX J21460478-0927054 


Abell 2400 


.0880 


II 


10/2000 


484 


1.05 


1.79 


2MASX J21574250-1 124402 


Abell 2401 


.0578 


II 


10/2000 


909 


1.07 


2.20 


2MASX J21582246-2006145 


Abell 2405 


.0366 




10/2000 


698 


1.13 


1.34 


2MASXJ21594215-1748019 ""^ 


Abell 2571 


.1081 


II 


11/1999 


570 


1.15 


1.46 


2MASX J23 183367-02 16356 


Abell 2721 


.1149 


II 


11/1999 


694 


1.01 


1.79 


2MASX J00055975-3443171 








10/2000 


808 


1.03 


1.33 




Abell 2730 


.120 


II 


11-12/1999 


818 


1.08 


2.38 


2MASX J00095644-3541015 = 


Abell 2804 


.1121 


I-II 


12/1999 


758 


1. 11 


1.98 


b 








10/2000 


758 


1.03 


1.48 




Abell 2811 


.1082 


I-II 


12/1999 


432 


1.30 


2.10 


2MASX J00420892-2832087 


Abell 2955 


.0945 


II 


10/2000 


695 


1.14 


1.41 


2MASX J01570191-1701234 


Abell 2969 


.1271 


I 


12/1999 


756 


1.03 


2.13 


2MASX J02033533-4106002 


Abell 2984 


.1044 


I 


10/2000 


870 


1.06 


1.57 


2MASX J021 12484-4017261 


Abell 3094 


.0683 


I-II 


10/2000 


745 


1.02 


1.60 


b 


Abell 3112 


.0759 


I 


12/1999 


530 


1.23 


1.87 


2MASX J03175766-4414175 


Abell 3166 


.1171 


I 


10/2000 


678 


1.04 


1.68 


2MASX J03464387-3248454 


Abell 3194 


.0927 


III 


10/2000 


878 


1.03 


1.80 


b 


Abell 3693 


.1237 




10/2000 


806 


1.09 


1.50 


2MASXJ20341991-3429387 "• 


Abell 3705 


.0883 


III 


10/2000 


813 


1.07 


1.41 


2MASX J20420438-3513067 


Abell 3727 


.1159 


III 


10/2000 


826 


1.04 


1.37 


2MASX J20593652-3629085 


Abell 3809 


.0616 


III 


10/2000 


790 


1.14 


1.37 


2MASX J21465904-4353564 


Abell 3920 


.1263 


I-II 


10/2000 


878 


1.23 


2.30 


2MASX J22492865-4053335 


Abell 4010 


.0963 


I-II 


11/1999 


709 


1.03 


1.73 


2MASXJ233 11412-3630405 


APMC 020 


.1105 




10/2000 


808 


1.05 


1.42 


2MASXJ00133506-3443086 "l 


Abell S0084 


.1087 


I 


10/2000 


763 


1.07 


1.25 


2MASX J00492282-293 1069 


Abell S0296 


.0699 


I 


12/1999 


514 


1.33 


2.30 


2MASX J02463700-42220 1 5 ^ 








10/2000 


899 


1.04 


1.44 




Abell S0952 


.0898 


II 


10/2000 


813 


1.16 


1.41 


b 


Abell SI 096 


.1049 


I-II 


10/2000 


760 


1.15 


1.49 


b 



^The BCG in this cluster has a complicated core with multiple nuclei, which are also visible in a Chandra image of the cluster. 

This BCG is not included in the statistical sample analyzed in this paper. 
'^ Abell 2405 is a superposition of two relatively poor systems at 11,000 and 27,000 km s~ . Here we focus upon the lower redshift group. 

No published Bautz-Morgan type. Abell 2405 and Abell 3693 each have two discrete redshift peaks. Published Bautz -Morgan types do exist for 
these clusters, but not for the redshift peaks that we are studying. 
^The listed value is the mean cluster redshift rather than the BCG redshift. 



The 1999 data for Abell S0296 are poor. The images contain a large gradient not present in any of the other data. 
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TABLE 2 
DE VAUCOULEURS PARAMETERS 



Cluster 


M 


CM 


re 


O-r,. 


h/a 


Cbja 


xl 


Abell 0122 


-24.98 


0.01 (0.02) 


29.9 


0.2(0.5) 


0.79 


0.01 


1.123 


Abell 1651 


-25.44 


0.01(0.02) 


57.7 


0.4(1.1) 


0.73 


0.01 


0.985 


Abell 2376 


-23.69 


0.01 (0.00) 


9.0 


0.1(0.1) 


0.88 


0.01 


0.987 


Abell 2400 


-24.77 


0.01(0.03) 


37.0 


0.8(2.9) 


0.60 


0.01 


0.997 


Abell 2401 


-24.28 


0.01 (0.00) 


11.4 


0.1(0.1) 


0.88 


0.01 


1.043 


Abell 2405 


-23.18 


0.01 (0.00) 


5.1 


0.1(0.0) 


0.85 


0.01 


1.002 


Abell 2571 


-24.66 


0.01(0.01) 


12.3 


0.1(0.2) 


0.94 


0.01 


1.046 


Abell 2721 (2000) 


-25.19 


0.01(0.01) 


40.7 


0.3(0.4) 


0.46 


0.01 


1.020 


Abell 2721 (1999) 


-25.18 


0.01 (0.03) 


39.1 


0.3(0.8) 


0.45 


0.01 


1.005 


Abell 2730 


-25.01 


0.01 (0.02) 


26.6 


0.2(0.7) 


0.86 


0.01 


1.066 


Abell 2811 


-24.78 


0.01(0.02) 


18.9 


0.2(0.4) 


0.68 


0.01 


1.022 


Abell 2955 


-24.54 


0.01(0.01) 


22.2 


0.2(0.4) 


0.77 


0.01 


1.068 


Abell 2969 


-24.92 


0.01(0.11) 


54.0 


0.7(10.1) 


0.71 


0.01 


1.022 


Abell 2984 


-24.65 


0.01(0.01) 


28.3 


0.2(0.4) 


0.77 


0.01 


1.109 


Abell 3112 


-25.47 


0.01(0.03) 


51.5 


0.3(1.2) 


0.52 


0.01 


0.987 


Abell 3166 


-24.05 


0.01(0.02) 


19.1 


0.2(0.4) 


0.80 


0.01 


1.014 


Abell 3693 


-24.15 


0.01 (0.00) 


9.1 


0.1(0.0) 


0.75 


0.01 


1.033 


Abell 3705 


-24.01 


0.01 (0.00) 


13.1 


0.1(0.0) 


0.86 


0.01 


1.041 


Abell 3727 


-23.79 


0.01 (0.00) 


12.8 


0.1(0.1) 


0.97 


0.01 


1.036 


Abell 3809 


-24.47 


0.01 (0.00) 


26.0 


0.1(0.2) 


0.80 


0.01 


1.006 


Abell 3920 


-24.79 


0.01(0.00) 


16.5 


0.1(0.0) 


0.78 


0.01 


0.994 


Abell 4010 


-25.00 


0.01(0.04) 


44.6 


0.4(2.1) 


0.68 


0.01 


1.019 


APMC 020 


-24.69 


0.01(0.00) 


22.4 


0.1(0.2) 


0.62 


0.01 


1.009 


Abell S0084 


-24.87 


0.01(0.00) 


30.2 


0.2(0.0) 


0.93 


0.01 


1.032 


Abell S0296 (2000) 


-24.97 


0.01(0.00) 


21.5 


0.1(0.1) 


0.69 


0.01 


1.053 


Abell S0296 (1999) 


-24.94 


0.01(0.01) 


21.6 


0.1(0.3) 


0.65 


0.01 


0.972 



Note. — The listed effective radius is Ve ^ yo^ and is in units kpc. The absolute magnitudes, 
which are total magnitudes that include extrapolation beyond r = 300 kpc, are not corrected for extinc- 
tion or k-dimming. Both corrections are typically <0.1 mag for galaxies in our sample. The column 
b/a gives the ratio of the minor to major axes. In the last column Xi, i^ the reduced x ■ For clus- 
ters observed during multiple years, the ensemble Xu value includes only the 2000 data. Parameter 
uncertainties corresponding to the Ic sky level uncertainties are given in parentheses. 
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TABLE 3 
Sersic Parameters 



Cluster 



M 



h/a <7i,i„ 



^x\ 



Abell 0122 


-26.13 


0.03(0.13) 


280.7 


Abell 1651 


-25.51 


0.01 (0.07) 


65.0 


Abell 2376 


-24.21 


0.02 (0.06) 


26.2 


Abell 2400 


-24.94 


0.02 (0.07) 


47.8 


Abell 2401 


-24.81 


0.01 (0.04) 


33.3 


Abell 2405 


-23.44 


O.OI (0.02) 


8.6 


Abell 2571 


-25.39 


0.02(0.16) 


58.8 


Abell 2721 (2000) 


-25.09 


0.01 (0.04) 


34.7 


Abell 2721 (1999) 


-25.12 


0.01 (0.07) 


35.5 


Abell 2730 


-26.88 


0.06 (0.39) 


1250.0 


Abell 2811 


-26.34 


0.06(0.31) 


522.0 


Abell 2955 


-26.01 


0.03(0.19) 


412.5 


Abell 2969 


-25.67 


0.03 (0.30) 


185.9 


Abell 2984 


-27.08 


0.06 (0.29) 


3019.2 


Abell 3112 


-25.95 


0.01(0.11) 


120.5 


Abell 3166 


-25.81 


0.08 (0.32) 


741.5 


Abell 3693 


-25.68 


0.10(0.25) 


489.7 


Abell 3705 


-24.64 


0.02 (0.05) 


44.4 


Abell 3727 


-25.16 


0.05 (0.20) 


207.1 


Abell 3809 


-24.53 


0.01 (0.02) 


29.3 


Abell 3920 


-25.14 


0.01 (0.04) 


33.0 


Abell 4010 


-25.89 


0.03 (0.25) 


223.4 


APMC 020 


-24.83 


0.01 (0.04) 


28.8 


Abell S0084 


-26.16 


0.04(0.13) 


315.8 


Abell S0296 (2000) 


-25.14 


0.01 (0.01) 


29.1 


Abell S0296 (1999) 


-25.04 


0.01 (0.04) 


25.9 



16.9 (:««-^) 
1.3 (!f:^) 

1.0 (+3-2) 

0.5 (!") 

0-1 di) 

2.4 (!?!;«) 

0.6 (!li) 

0.7 (!^|) 

177 4 (+1759.3) 

75.4 (f if) 



I 3 (+I42.8-, 



2 5 (+2**'") 
143.8 (||^) 
142.3 ( J|2) 



0.3 C 



■+\1\ 



0.5 (+[«) 

24 1 (+'**-5) 
■^^■^ '•-60.2-' 

0.2 (+0|) 
0.3 Cjl) 



9.60 
4.26 
7.47 
5.02 
7.52 
5.78 
9.24 
3.65 
3.11 
14.36 
13.65 
11.52 
111 
14.81 
6.08 
14.16 
17.74 
7.40 
11.75 
4.29 
6.38 
1.16 
4.67 
9.53 
4.70 
4.54 



0.14 (; 



-0.16 
+0.31 



0.11 ( 

0-06 csj:?^ 

0.04 (;0;| 
03 r+''' ' 
13 r+' '■♦ 

^-0.18 

37 (+2*9 

MSI 

0.41 Cfl 

0-17 (S 
0.12(:i;«« 

0.28 C'l 

0.04 (!; 

52 (+'''^^ 
M.5Q 

1.03 («]= 

0.10 (!»i 

0.32(:>;15 
02 (+''*» 
0.08 (!«^i) 
0.10 (t?|) 
0.04 0.\\) 
0.17 (-™-56) 
0.01 (f |) 

0.03 (!«:{5) 



0.73 
0.73 
0.86 
0.61 
0.89 
0.84 
0.95 
0.46 
0.45 
0.84 
0.65 
0.73 
0.70 
0.72 
0.50 
0.77 
0.80 
0.85 
0.93 
0.80 
0.77 
0.65 
0.62 
0.93 
0.68 
0.65 



0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 



-6032 
-38 
-2073 
-195 
-14411 
-5774 
-4580 
-84 
-19 
-5253 
-3351 
-9107 
-390 
-10569 
-3406 
-2693 
-2535 
-3419 
-2494 
-165 
-1296 
-3198 
-270 
-3605 
-3338 
-302 



Average 



-3680 



Note. — In the last column Axd = Xwrsic ~ XdeVaucouleur^- F°i" clusters observed during multiple years, the average and ensemble values 
are computed using only the 2000 data. Parameter uncertainties corresponding to the Icr sky level uncertainties are given in parentheses. 
As in Tablelzl the absolute magnitudes are not corrected for extinction or k-dimming. These absolute magnitudes are total magnitudes that 
include extrapolation beyond r = 300 kpc. 
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TABLE 4 
2-deV Parameters 



Cluster 



M,a,al 



|A6»| 



o"|Ae| 



M 



Te 



bja a^i„ 



^x\ 



Ax? 



Abell 0122 
Abell 1651 
Abell 2376 
Abell 2400 
Abell 2401 
Abell 2405 
Abell 2571 
Abell 2721 (2000) 
Abell 2721 (1999) 
Abell 2730 
Abell 2811 
Abell 2955 
Abell 2969 
Abell 2984 
Abell 3112 
Abell 3166 
Abell 3693 
Abell 3705 
Abell 3727 
Abell 3809 
Abell 3920 
Abell 4010 
APMC 020 
Abell S0084 
Abell S0296 (2000) 
Abell 80296(1999) 



Average 
Average/dof 



-25.66 
-25.84 
-23.97 
-25.17 
-24.62 
-23.91 
-25.42 
-25.20 
-25.19 
-26.14 
-25.70 
-25.31 
-25.58 
-25.87 
-25.82 
-24.73 
-25.63 
-24.38 
-25.14 
-24.76 
-25.08 
-25.71 
-25.12 
-25.38 
-25.20 
-25.05 



/+0.03 
'■-0.04 



/+0.I8 
'.-0.32 

/+0.03 
'-0.03 

,-+0.17 
'•-0.28 



/+0.15 
'.-0.lt) 

/+0.02 
'.-0.02 

/+0.03 
'--0.04 



/+0.24 
'--0.30 



/+0.22 
'--0.22 



/+0.08 
'.-0.08 

/+0.26 
'.-0.44 

,-+0.01 
'.-0.09 



,-+0.07 
'-0.07 



/+0.04 
'--0.08 



.+0.2K 
^-0.28' 



/+U.0 3 % 
'--n oi ' 



'.-0 04' 



76.5 1.1(0.5) 

9.2 1.3(0.5) 

68.5 13.2(5.5) 
10.0 1.5(2.5) 

24.4 3.2 (0.4) 

43.2 2.5 (0.3) 

83.6 0.8 (0.5) 

1.1 1.3(0.4) 
9.7 1.4(0.5) 

40.3 5.0 (4.6) 
8.7 3.3 (4.9) 

70.0 1.9(2.2) 

55.5 9.4 (28.9) 
64.5 8.0 (5.0) 

42.3 15.5 (29.5) 

82.4 2.1(0.1) 

71.5 5.4(0.6) 
2.5 4.1 (2.6) 

88.3 7.7(1.5) 

4.2 1.4(0.6) 
0.5 3.3 (0.4) 
9.4 3.2 (7.4) 
2.4 1.0(0.3) 

73.6 3.0(1.3) 

7.1 0.8(0.9) 

24.2 3.2 (4.5) 



-25.56 
-22.97 
-24.95 
-25.21 
-23.87 
-21.33 
-24.61 
-24.19 
-24.48 
-22.30 
-23.38 
-22.89 
-25.11 
-23.89 
-25.02 
-23.14 
-25.04 
-22.97 
-26.00 
-23.85 
-25.53 
-23.59 
-25.22 
-22.61 
-25.35 
-23.80 
-25.79 
-23.04 
-25.74 
-22.95 
-24.65 
-21.92 
-25.46 
-23.56 
-24.27 
-21.85 
-24.97 
-23.07 
-23.85 
-24.14 
-24.80 
-23.48 
-25.56 
-23.52 
-24.34 
-24.39 
-25.33 
-21.94 
-24.89 
-23.67 
-24.82 
-23.25 



-0.05 
+0.48 



-M 



-0.11 
+0.02 -, 
-0.02.' 
+0.10-, 

%\ 



m 



0.01 
0.03 
0.07 
0.01 
0.01 
0.08 
0.04 
0.08 
0.01 
0.03 
0.02 
0.01 
0.01 
0.02 
0.03 
0.14 
0.04 
0.21 
0.03 
0.03 
0.03 
0.03 
0.01 
0.03 
0.03 
0.08 
0.01 
0.02 
0.01 
0.04 
0.02 
0.04 
0.05 
0.02 
0.01 
0.07 
0.05 
0.03 
0.03 
0.03 

0.0 

0.10 (! 

0.02 (til) 



-0.32.' 
+0.19-, 
-0.16^ 
+0.22 -, 
-0.23.' 
+0.22 -, 

-0.08^ 
+0.09 -, 

+o:i? 

-0.49^ 
+0.21 -, 

-8:8? 

-0.04^ 
+0.06 -, 
-0.06' 
+0.16-1 

+o:8L 

-0.08' 
+0.20-, 

;8:?i\ 
;o°o§t 

-o.oo' 

+0.18-, 



mi 

-0.16' 
+0.02 -, 

-ml 

'--Q.Q^'' 



0.06 (+« " 
0.04 ( 



0.03 (! 
0.01 (t_ 



0.04 (!<!«|j) 

'--Q.0Q-' 

0.05 (!«• j«) 
0.04 (^p 
0.22 (^2?) 



107.9 
4.7 
477.2 
48.2 
18.4 
1.0 
184.6 
20.4 
27.6 
1.8 
93.9 
3.8 
107.7 
5.4 
47.9 
18.1 
45.8 
16.1 
245.9 
8.1 
139.5 
5.5 
93.5 
3.4 
232.0 
19.7 
202.7 
6.1 
102.1 
8.3 
69.2 
1.7 
314.5 
4.9 
30.8 
2.1 
215.2 
6.1 
126.8 
20.8 
45.6 
5.4 
186.3 
13.7 
180.2 
17.6 
72.0 
2.0 
46.1 
10.0 
33.9 
8.6 



2.9 


(.+9.3 X 


1 


f+U.2-, 


53.1 ( 
0.5 


r+8^gf 

Lo.2^ 


0.6 


f+2.0-, 
Ll.o' 


1 


(-+0.2-I 


44.9 ( 


Ml. 
75 2 ^ 


1.0 


f+faSx 
L1.4' 


0.4 


f+2.3-| 


1 


f+U.2-, 


4 7 ( 


;1^, 


0.1 

5 Of 




1 


r+^'.r 


0.8 




1.6 


(-+1.2-, 


0.9 
1 9 




13.5 ( 


+/219-, 

-78.4 > 


0.2 


f+1.2-. 


8.6 ( 


+59.7-, 
-44.4' 


0.2 


r+0.8-, 

Li.o' 


2.3 ( 


+11.3-, 
-9.8 ' 


1 


f+0.3-. 


26.7 ( 


+1'o^;6. 

-83.2 > 


0.8 


Llj' 


5.4 ( 
0.1 


+23.2-, 
Lo.o' 


1.7 ( 


+12.4-, 
-9.4 ' 


0.2 
-76 


f+0.9-, 

A 


0.1 


Lo.l' 


74 4 


f+bl.5-, 


1 ( 


^^,\ 


1.0 ( 


^t^9 
-2.6 ' 


0.1 


f+1.5-. 


17.1 


/+54.6-, 
'-^3.7' 


1 


(-+0.2-, 


17.2 


Ma. 

'-22.0' 


0.4 

3.2 


r+0. 1 -, 

■di) 


0.^ 


Kli) 


10.7 


,-+107.4 
'--58.1 


O.f 


ir+2-8^ 

'+i5§^l 

'-7L3 


34.4 


0..i 


1 ,+0.9 -. 


l.f 


0. 


/+0.2^ 
'--0.2-' 


0.: 


/+7.5^ 

A 

'--0 3'' 


2.1 




0.^ 


' ^^-0.6' 



0.53 


0.01 


-9189 


-3156 


0.63 


0.01 






0.27 


0.01 


-328 


-289 


0.78 


0.01 






0.83 


0.01 


-2197 


-123 


0.89 


0.04 






0.41 


0.02 


-339 


-144 


0.67 


0.01 






0.88 


0.01 


-15350 


-939 


0.87 


0.01 






0.81 


0.01 


-7687 


-1912 


0.83 


0.01 






0.52 


0.01 


-7247 


-2667 


0.68 


0.01 






0.43 


0.01 


-810 


-725 


0.44 


0.01 






0.43 


0.01 


-391 


-372 


0.42 


0.02 






0.72 


0.01 


-5382 


-129 


0.92 


0.01 






0.49 


0.01 


-3556 


-204 


0.85 


0.01 






0.57 


0.01 


-10927 


-1820 


0.79 


0.01 






0.39 


0.01 


-841 


-450 


0.93 


0.01 






0.52 


0.01 


-11873 


-1303 


0.95 


0.01 






0.41 


0.01 


-5666 


-2260 


0.96 


0.01 






0.61 


0.01 


-3355 


-662 


0.52 


0.03 






0.84 


0.02 


-3713 


-1178 


0.68 


0.01 






0.93 


0.01 


-3522 


-103 


0.62 


0.02 






0.68 


0.02 


-2953 


-458 


0.95 


0.01 






0.33 


0.02 


-1096 


-931 


0.91 


0.01 






0.85 


0.01 


-1421 


-124 


0.68 


0.01 






0.47 


0.01 


-3579 


-381 


0.88 


0.01 






0.32 


0.02 


-737 


-467 


0.67 


0.01 






0.85 


0.01 


-4027 


-422 


0.66 


0.03 






0.57 


0.01 


-7429 


-4091 


0.78 


0.01 






0.59 


0.01 


-693 


-390 


0.69 


0.01 










-4718 


-1039 






-786 


-208 



Note. — In the last two columns Ax^ = Xi-deV ~ Xde Vaucouieurs ^"'^ ^X^ = Xi-deV ~ Xsersic- Parameter uncertainties corresponding to the 1 a sky level uncertainties 
are given in parentheses. As in Tablesl^andj^ the absolute magnitudes are not corrected for extinction or k-dimming. These absolute magnitudes are total magnitudes, 
which include extrapolation beyond r = 300 kpc. 
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TABLE 5 
Abell 1651 Sersic Parameters 



Gonzalez et al. (2000) 59.6±1.1 4.3±0.2 
This work 65.0^*^ 4.26^^ [^ 





TABLE 6 

Average Absolute Magnitudes 








Model 


MtoicI f^MT^^i 


^300 o-M,oo 


M50 


°"«50 


Mio 


O-M,0 


de Vaucouleurs 

Sersic 

2-deV 


-24.60 0.56 
-25.43 0.83 
-25.23 0.59 


-24.56 0.53 
-24.93 0.50 
-24.96 0.47 


-24.20 
-24.26 
-24.26 


0.39 
0.36 
0.35 


-23.14 
-23.11 
-23.11 


0.30 
0.28 
0.29 



Note. — The average absolute magnitudes and the dispersion in these magnitudes for our BCG sample. 
The column Mjotai conesponds to the total absolute magnitudes hsted in Tables|2|-|4| The columns M300, 
M50, and M\Q correspond to the model magnitudes within circular metric apertures of radii 300, 50, and 10 
kpc, respectively. The above table does not include extinction or k-corrections, but both effects are small. 
Accounting for these factors decreases the dispersions by only 0.01 mag. The impact on the mean absolute 
magnitude can be approximated by assuming a mean extinction of 0.05 mag and a mean k-con"ection of 
0.07 mag. 



TABLE 7 

Impact of PSF on Derived 2-deV Profile for Abell 0122 



PSF 


FWHM(") 


Md 


t-e.i 


(blah 


A»e.2 


rj,2 (kpc) 


{b/a)2 


|A6»| 


Ax^ 


Abell 0122 


1.38 


18.77 


5.9 


0.63 


19.67 


147.8 


0.53 


76.5 


_ 


Abell 2376 


1.33 


18.83 


7.4 


0.67 


19.74 


157.9 


0.52 


72.6 


51.7 


Abell 3727 


1.37 


18.81 


6.6 


0.63 


19.71 


152.6 


0.52 


75.3 


82.4 


Abell SI 096 


1.49 


18.80 


6.6 


0.63 


19.71 


153.8 


0.53 


73.3 


-23.9 



Note. — Ax is the difference between the x using this PSF and the x using the conect PSF for the Abell 
0122 field. 



